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Abstract 

In this paper we propose a new version of differential transform method (we 
shall call this method as a-parameterized differential transform method), 
which differs from the traditional differential transform method in calculat¬ 
ing coefficients of Taylor polynomials. Numerical examples are presented 
to illustrate the efficiency and reliability of own method. The result reveal 
that a-Parameterized differential transform method is a simple and effective 
numerical algorithm. 
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1. Introduction 

Many problems in mathematical physics, theoretical physics and chem¬ 
ical physics are modelled by the so-called initial value and boundary value 
problems in the second-order ordinary differential equations. In most cases, 
these problems may be too complicated to solve analytically. Alternatively, 
the numerical methods can provide approximate solutions rather than the 
analytic solutions of problems. There are various approximation methods 
for solving a system of differential equations, e.g. Adomian decomposition 
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method (ADM), Galerkin method, rationalized Haar functions method, ho- 
motopy perturbation method (HPM), variational iteration method (VIM) 
and the differential transform method (DTM). 

The DTM is one of the numerical methods which enables to hnd ap¬ 
proximate solution in case of linear and non-linear systems of differential 
equations. The main advantage of this method is that it can be applied di¬ 
rectly to nonlinear ODEs without requiring linearization. The well known 
advantage of DTM is its simplicity and accuracy in calculations and also wide 
range of applications. Another important advantage is that this method is 
capable of greatly reducing the size of computational work while still accu¬ 
rately providing the series solution with fast convergence rate. With this 
method, it is possible to obtain highly accurate results or exact solutions for 
differential equations. The concept of the differential transform method was 
hrst proposed by [l^, who solved linear and nonlinear initial value prob¬ 
lems in electric circuit analysis. Afterwards, Chiou and Tzeng applied 
the Taylor transform to solve nonlinear vibration problems, Chen and Ho 
developed this method to various linear and nonlinear problems such as 
two point boundary value problems and Ayaz applied it to the system 
of differential equations. Abbasov et al. [l| used the method of differential 
transform to obtain approximate solutions of the linear and non-linear equa¬ 
tions related to engineering problems and observed that the numerical results 
are in good agreement with the analytical solutions. In recent years many 
authors has been used this method for solving various types of equations. For 
example, this method has been used for differential-algebraic equations 14, 
partial differential equations ,0,0i0, fractional differential equations [l^ 
and difference equations [l^. In jol. Iiol Ed| . this method has been utilized 
for Telegraph, Kuramoto-Sivashinsky and Kawahara equations. Shahmorad 
et al. developed DTM to fractional-order integro-differential equations with 
nonlocal boundary conditions [l2j | and class of two dimensional Volterra in¬ 


tegral equations [13|- Abazari et. al. applied this method for Schrodinger 
equations 141. Different applications of DTM can be found in 0, @] . Even 
if the differential transform method (DTM) is an effective numerical method 
for solving many initial value problem, there are also some disadvantages, 
since this method is designed for problems that have analytic solutions (i.e. 
solutions that can be expanded in Taylor series). 

In this paper we suggest a new version of DTM which we shall called 
a-parameterized differential transform method (a-p DTM) to solve initial 
value and boundary value problems, particularly eigenvalue problems. Note 
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that, in the special cases a = 1 and a = 0 of the a-p DTM reduces to the 
standard DTM. 


2. The classical DTM 


In this section, we describe the dehnition and some basic properties of 
the classical differential transform method. An arbitrary analytic function 
f(x) can be expanded in Taylor series about a point x = Xq as 


/(^) = 


X 


XoY 


k r Jk 


k=0 


k\ 


dx^ 


X=Xq 


( 1 ) 


The classical differential transformation of f(x) is dehned as 

( 2 ) 

X=Xq 

Then the inverse differential transform is 




d^f{x) 

dx^ 


fi^) = - Xo)^F{k). 


(3) 


k=0 


The fundamental mathematical operations performed by differential trans¬ 
form method are listed in following 

i) If /(x) = g{x) ± h{x) then F{k) = G{k) ± F[{k) 

ii) If /(x) = ag{x), a G M, then F{k) = aG{k) 

hi) If /(x) = then F{k) = {k + 1)(A: -|- 2)...(A: -|- m)G{k + m) 

1 for k = m 
0 for k m 


iv) If /(x) = x”^ then F{k) = 5{k — m) = 


v) If /(x) = g{x)h{x) then F{k) = ^ H{m)G{k 

m=0 


m] 


3. a-Parameterized differential transform method (a-p DTM) 

In this section we suggest a new version of classical differential transform 
method by following. 

Let J = [a, 6] C M be an arbitrary real interval, /:/—)■ R is an inhnitely 
differentiable function (in real applications it is enough to required that /(x) 


3 



is sufficiently large order differentiable function) , a G [0,1] any real pa¬ 
rameter and N any integer (large enough). Let us introduce the following 
notations 


^a(/;fc) 


k\ 


Db{f;k) 


k\ 


(4) 


D(/, a; k) := aD^', k) + {I - a)D,{f. k), 


a e [0,1], k eN 


(5) 


Definition 1. The sequence 

(D„(/)) :=(D(/,a;l),D(/,a;2),...) 

is called the a-P transformation of the original function f(x). The so-called 
,,differential inverse” transform of Da{f) we dehne as 

CO 

Eo.{DM)) ■.= J2D{f,a-k){x-x^)'^ (6) 

k=0 

if the series is convergent, where Xa = cta + {1 — a)b. 

The function fa{x) dehned by equality 

Ux) := E^{DM)) 

we called the a-parameterized approximation of the function f(x). 

Remark 1. In the cases of a = 1 and a = 0 the a-p differential trans¬ 
form ^ reduces to the classical differential transform ^ at the points x = 
and x = h respectively. Namely for a = 0 and a = 1 the equality 
fa{x) = f{x) is hold. 

Remark 2. For practical application, instead of fa{x) it is convenient to in¬ 
troduced N-term a-parameterized approximation of the function fa{x) which 
we shall define as 


N 

Ia,N{x) := Ea,N{Do,{f)) := ^ D{f, a; k){x - x^)^ 

k=0 


(7) 
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Theorem 1. If f(x) is constant function then fa{x) = f{x) and fa,Nix) = 
fix) for each N. 

Proof. The proof is immediate from Definition [1] and Remark [2l 

Theorem 2. If fix) = cgix), c G R, then Daif) = cDaig) and /q(x) = 
cgaix). 

Proof. By applying the well-known properties of classical DTM we get 

Dif,a;k) = aDaif;k) + il-a)Dbif;k) 

= acDaig;k) + il - a)cDbig;k) 

= ciaDaig;k) + il - a)Dbig;k)) 

= cDig,a]k) ( 8 ) 

Consequently Daif) = cDaig), from which immediately follows that /q,(x) = 
cgaix). 

Theorem 3. If fix) = gix)±hix) then Daif) = Daig)IcDaih) and faix) = 

cfaix) ± haix). 

Proof. By using the dehnition of transform ([5]) 

Dif,a-k) = aDaif-k)±il-a)Dbif-,k) 

= aDaig + h-,k) ±il - a)Dbig + h-,k) 

= Dig,a;k) ± Dih,a;k) (9) 

Consequently Daif) = Daig) hDaih), from which immediately follows that 
faix) = gaix) ± haix). 

Theorem 4. Let fix) = and m G N. Then 

Dif"^\ a; k) = a]k + m) 

and 

7M(x) = ^ UlL^-p^Dif, a;k + m)ix- Xa)^ 

k=0 

where = aa -|- (1 — Q)h. 
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Proof. We have from definition ([S]) 

= a{k + l){k + 2)...{k + m)Da{f; k + m) 

+ (1 - q;)(A: + l)(fc + 2)...{k + m)Db{f; k + m) 

= {k + l){k + 2)...{k + m){aDa{f; k) + {1 - a)Db{f; k + m)) 

= ^^^^jp^D{f,a;k + rn) 

Thus we get a; k) = (^+^)' £)(y^ q,j k+m). Using this we hnd fjr'\x) = 

OO 

a;k + m){x - x^Y- 

k=0 

Theorem 5. Let f{x) = x^, m G N. Then 

^ ^ + (1 ~ a)b"^~Y for k < m 

1 for k = m 
0 for k > m 

Proof. Let k < m. By using the dehnition of the transform ([5]) we have 
DU,a-k) = aDYf;k) + {l-a)Db{f-,k) 

= + (1 

= ^ ^ + (1 - ( 10 ) 

The equalities D{x‘^,a]m) = 1 and D{x‘^,a]m + s) = 0 for s > 1 is 
obvious. 

k 

Theorem 6. If f{x) = g{x)h{x) then D{f,a]k) = [(^IIa{g]' 02 )Da{h] k — 

m=0 

m) + (1 — a)Db{g] m)Db{h] k — m)] 

Proof. By using the dehnition of transform given in Eq. (jS]) we have 
D{f,a-,k) = aDa{f-,k)+ m{l-a)Db{f-,k) 

= aDa{gh;k) + {1 - a)Db{gh;k) 

= wX(f) E (™ 

m=0 I m=0 I 

k 

= '^[aDa{g]m)Da{h]k-m) + {l- a)Db{g]m)Db{h]k-m)]. (11) 

m=0 
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4. Justification of the a-p DTM 

In order to show the effectiveness of o-p DTM for solving boundary value 
problems , examples is demonstrated. 

Example 1. (Application to boundary-value problem) Let us consider the 
differential equation 

£y := y''{x) + y^^y{x) = 0, a: G [0,1], // G M. (12) 

with the boundary conditions 


1 /( 0 ) = 0 , 1 /( 1 ) = 1 . 

Exact solution for this problem is 

sin /ix 


y{x) = 


sm fi 


(13) 


(14) 


Applying the A^-term a-p differential transform of both sides (IT^ and (IT^ 
we obtain the following a-parameterized boundary value problem as 

{^y)a,N = 0a,N, 2/o,Af(0) = 0a,N, 2/o;Af(l) = la,Ar- (15) 

By using the fundamental operations of a-p DTM we have 

(16) 

^ {k + l){k + 2) ^ ^ 

The boundary conditions given in flljp can be transformed as follows 

N N 


y< 


,,n{ 0) = '^D{y,a]k){a-l)’^ = 0 and i/a,jv(l) = D{y, k)a’^ = 1(17) 


k=0 


k=0 


Using (1T6|) and (ITT)) and by taking N = 5, the following a-p approximate 
solution is obtained 


ya{x) = A+{x-Xa)B- 

y^{x — XaYB 


y?{x — XaYA y?{x — Xa^B ^ jj'^x — XaYA 


6 


24 


120 


+ 0{x^) 


(18) 
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where Xq, = (1 — a), according to ([S]), D{y,a-,0) = A and D{y,a-,1) = B. 
The constants A and B evaluated from equations in fll6l) as follow 


A = 


X 


(2880 (xq/U® — + Oxg/i® + — Ax^y^ — 480/i^ + Xq/t® + 2880 


12xo/u® + 40xq//® — 60xo/i® + ASxlf/ — lOxg/i®) 


{xq 


y'^xl ^ y^xl ^ 


6 


120 


(19) 


and 


B = 


X 


(2880 (x^/ - 4x^/ + 6x®/i® + 24/i^ - 4x^/i® - 480yU^ + x®/ + 2880 


12xo/i° + 40xQ/i° — OOxq/x” + 48xo/i“ — IOxq//®) ) 


2 2 4 4 

^ 2 24 ' 


( 20 ) 


Example 2. (Application to eigenvalue problems)We consider the following 
eigenvalue problem 


y" + Xy = 0, X e [0, l] 

4 ii 2 /( 0 ) + Ai2y'{Q) = 0 
421 |/( 1 ) + 4221 /^( 1 ) = 0 

Taking the a-p differential transform of both sides fl2ip we find 

D{y" + Ay, a; A;) = (A; + 1)(A; + 2)D(y, a;k + 2) + ^D(y, a; k) = 0. 

Thus the following recurrence relation is obtained 

r,/ 7 , o^ _ ^D{y,a-k) 

(A: + l)(A: + 2)' 

Using definition of a-p differential transform we get 


( 21 ) 

( 22 ) 

(23) 


(24) 


(25) 


y„(x) = ^ D(i/, a; k)(x - x^f 

k=0 

oo 

y'c(x) = ^ kD{y, a; k){x - 


(26) 


fc =0 


(27) 







Consuquently 


5„(0) = D(y, a; k){a-l)'‘ = 5^(-l)‘C(!,, a; fc)(l - a)‘ (28) 

k=0 k=0 

OO CXD 

y'ai^) = a- k){l - af-^ (29) 

fc =0 fc =0 

Thus the boundary condition given in fl2^ can be transformed as follows 

CXD 

^ii2/«(0) + 24i2y'„(0) = ^^(a4ii(Q; — 1)^ + kAi2{a — 1)^ ^)D{y, a; k) = 0(30) 

k=0 

Similarly we have 


k=0 

and 

CXD 

y>^{l) = Y,kD{y,a;k)a’^-^ (32) 

k=0 

In this case the boundary condition given in 0231) can be written as follows 


^21?/a(l) + ^222 /'q-( 1) — ^^(^2lCK^ + D {V, Oi] k) — 0 (33) 

fc =0 


Let D{y, a;0) = A and D{y, a; 1) = B. Substituting these values in fl25l) 
we have the following recursive procedure 


D{y,a;k) 


( 2 ^!) ’ 
B{-\V 

(2£+l)! ’ 


for k = 2i 
for k = 2i + 1 


Substituting (IM)) in (l30|) and (l33|) we hnd 


(34) 


^{^(24ii(a; - 1)^^ + 2Mi2(a - 1)^^ ^) J, } 
e=o ^ '' 

°° (—\V 

+ B{J2{An{a - 1)^^+^ + (2£ + l)^2(a - = 0 (35) 

e=o I + J- 
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and 


e=o ^ '' 

°° (—\Y 

+ B{^{A2ia'^^~^^ + {2i + 1 )^ 220 ^^) ^2i-{- 1 ) 1 ^ ^ 

respectively. In this case we have a linear system of the eqnations with 
respect to the variables A and B as 


APn(A) + 5Pi2(A) = 0 (37) 

v4i^2i(A) + BP22{X) = 0 (38) 

where Pii(A) := + 2 Mi2(q; - 1)^^"^)^^, ^i 2 (A) := 

~ 1)2^+^ + (2£+ 1)Ai2(q! - l)^'^) (2^+1)!, -P2 i(A) := + 

2M22a^'^“^)^^ andP22(A) := + (2f +1)^22^^^ (27+1)! • 

the system flTr)) - fl38|) has a nontrivial solntion for A and B the characteristic 
determinant is zero i.e. 


P(A) 


Pii(A) Pi2(A) 
P2i(A) P22(A) 


0 . 


Thns we have characteristic eqnation for eigenvalnes by a-p DTM Now, let ns 
hnd the exact eigenvalnes and eigenfnnctions of the Stnrm-Lionville problem 
fl^ - fl23]) . The general solntion of eqnation fl2^ have the form 


y{x) = C cos fix + D sin fix (39) 

where A = and (7, D are arbitrary constants. Applying the bonndary 
conditions ([21]), (122]) we get 


AiiC + fiAi2D — 0 (40) 

(A 21 cos —fiA22 sin fi)C + (A 21 cos —fiA22 sin fi)D = 0 (41) 

Becanse we cannot have (7 = ^ = 0, this implies 

(A 11 A 21 + fX^A 12 A 22 ) sin/i — fi{^Ai 2 A 2 i — A 11 A 22 ) cos/i = 0. (42) 
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This is transcendental eqnation which is solved graphically. Let /i = G 
N are points of intersection of the graphs of the fnnctions 


y = (^ 11^21 + h^^i 2 ^ 22 ) sin/i and (43) 

y — h(^ 12^21 ~ ^ 11 ^ 22 ) cos fi (44) 

The eigenvalnes and corresponding eigenfnnctions are therefore given by 

A„ = and yn{x) = Cn cos sin n G N. (45) 

Now we consider special case of the Stnrm-Lionville problem fl2^ - fl2l|) . 

y" + Xy = 0 (46) 

y{0) + y\0) = 0 (47) 

y{l)-y'{l) = 0. (48) 


The eigenvalnes of this problem are determined by the eqnation 

2yL 

tan^u = -- 

1 — 

This eqnation can be solved graphically by the points of intersections of the 
graphs of fnnctions 



1 

y = tan/i and y= - - 

1 — 


(50) 


as shown by the seqnence (//„) in Fignre 1. 

The eigenvalnes of the considered problem are given by A„ = and 
corresponding eigenfnnctions are given by 


yn{,x) = Cn COS + Dn sin y,nX, n G N. 


(51) 


Taking the a-p differential transform of both sides the eqnation fl4^ the 
following recurrence relation is obtained 


D{y,a;k + 2) = 


XD{y,a;k) 
{k + 1) (A) + 2) 


(52) 
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Figure 1: 


Applying N-term a-p differential transform the boundary conditions (l47|l - 
(H8|l are transformed as follows: 

N 

ya(0) + y'aW = a-,k) = 0 (53) 

k=0 

N 

ya(l) - 2/'a(l) = - ka^~^)D{y, a]k) = 0 (54) 

A:=0 


By using fl5^ . fl5^ and fl54)) we obtain the following equalities (for N=6) 
A 


1 + ((a - 1)^ + 2 (q! - 1)) ^ 21 ^ + ((“ “ 1)^ + 4(tt - 


+ ((« — 1)® + 6(a: — 1)^) 


-A3) 


6 ! 

A^i 


+ B 


a + ((a — 1)3 + 3(a — 1)^) 


2^(-A) 


3! 


+ ((«-l)3 + 5(a-l)^)- 


= 0 


(55) 


and 


A 

B 


1 + (a^ - 2a) ^ + {a^ - da^)^ + (a® - ba®)^—^ 

21 4! o! 

(a — 1) + (q;3 — 3q;^)-^—^ + (a® — 50;"^)^ = 0 

3! 5! 


(56) 
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Since the system fl55p - fl56p has a nontrivial solution for A and B the charac¬ 
teristic determinant is zero i.e. 


a(A) 


aii(A) ai2(A) 

021(A) a 22 (A) 


where an = 1 + ((a - 1)^ + 2(a - l))^p + ((a - 1)"^ + 4(a; - 1)^)^ + ((a - 

1)® + 6(q! - ai 2 = a + ((a - 1)^ + 3(a - 

+ ((a - 1)® + 5( q ! - 021 = 1 + - 2 a)^p + (a‘^ - + (a® - 

6 q ;^) ^ , 022 = (a — 1) + (a^ — + (a^ — 

Taking a = | we have the following algebraic equation for approximate 
eigenvalues: 

A 11A2 89A3 299A^ llA'^ 

“ 6 120 ~ 15360 ^ 2211840 ~ 9830400 “ ^ 

This equation can be solved by various numerical methods. 


5. Analysis of the method 

In this study we introduce a new version of classical DTM that will extend 
the application of the method to spectral analysis of boundary-value problems 
involving eigenvalue parameter, which arise from problems of mathematical 
physics. Numerical results reveal that the a-p DTM is a powerful tool for 
solving many initial value and boundary value problems. It is concluded that 
comparing with the standard DTM, the a-p DTM reduces computational cost 
in obtaining approximated solutions. This method unlike most numerical 
techniques provides a closed-form solution. It may be concluded that a-p 
DTM is very powerful and efficient in hnding approximate for wide classes 
of boundary value problems. The main advantage of the method is the fact 
that it provides its user with an analytical approximation, in many cases 
an exact solution, in a rapidly convergent sequence with elegantly computed 
terms. 
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